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SUMMARY 



Following the development of the well-known differen- 
tial equations of the problem and their resolution for 
failure in tension, the bending (transverse) oscillations 
of an originally not quite straight bar hinged at both 
ends and subjected to a constant longitudinal force (shock 
load) are analyzed. To this end the course of the bar form 
is expanded in a sinusoidal series, after which the inves- 
tigation is carried through separately for the fundamental 
oscillation and the (n-l)th higher oscillations. 

The analysis of the fundamental oscillation distin- 
guishes three cases: shock load lower, equal to, or higher 
than the Eulerian load. 

The investigation of the (n-l)th higher oscillation 
also distinguishes between shock load smaller, equal to, 
or greater than the (n-l)th stability limit, although on- 
ly the first case is of practical significance. 

Shock loads in buckling are divided into the period 
of actual shock and the period of free oscillations follow- 
ing the actual shock. 

The investigation leads to functions which are propor- 
tional to the maximum stresses in time and space due to 
the shock stresses in buckling. These functions are then 
compared for the case of shock load lower than Eulerian 
load with the maximum stresses in static load. It is found 



* " Stossart ige Knickbeanspruchung schlankor Stabe im elas- 
tischen Bereich bei beiderseits gelenkiger Lagerung." 
Luftf ahrtf orschung, July 6, 1933, pp. 55-64. 
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that the former are smaller for short shock periods and 
vice versa; that is, in the extreme case, twice as high as 
the latter. 

From a comparison of the functions decisive for the 
maximum stresses, it appears that the Eulerian load may "be 
safely exceeded in shock-like buckling stresses, provided 
the shock period is sufficiently short; further, that, 
whereas the stresses under shock load above the Euler load 
show an unrestricted increase with the shock period, the 
stresses in shock loads below the Euler load roach an up- 
per limit which is not exceeded during any shock period. 

The report closes with an analysis of the interdepend- 
ence between the shock stress in buckling and the shock 
impulso / P dt. It is found that, contrary to common be- 
lief, the stress with equal shock impulso is sensibly af- 
fected by the shock period. For that reason the determi- 
nation of the stress stipulates not only the time integral 
/ P dt, but also the shock force and the shock period - 
a fact which is of essential importance from the experi- 
mental point of view. 

I. INTRODUCTION 



The analysis of static buckling stresses affords, as 
is known, a problem in stability. It poses and answers 
the question up to what limit the compression may be in- 
creased for given bar dimensions without exceeding the 
range within which an unequivocally definable condition of 
equilibrium exists. Several equilibrium conditions are 
possible after this boundary has been exceeded. On ap- 
proaching the stability limit the rise of the deformation 
is such that the bar usually loses its carrying capacity 
before reaching the equilibrium condition. For this rea- 
son, the determination of the stability limit is of de- 
cisive importance. 

Contrariwise, stressing a bar suddenly in buckling, 
the suddenness being the short-time interval between load 
change and loading period, as shown in this report, the 
stability limit is no longer as significant as in the stat- 
ic case, and may be safely exceeded, provided the shock 
period is so short as to leave the bar no time to deform 
as would correspond to the static equilibrium condition. 
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From this it follows that the calculation of the def- 
ormations and stresses with respect to time is the primary 
issue rather than the determination of the stability lim- 
it when analyzing the shock load in buckling. 



II, NOTATION 



L 


kg, 


force component parallel to bar axis (x) . 


D 


kg, 


force component at right angles to bar axis 

(y). 


X=- P 


kg, 


force component parallel to axis x. 


Y 


kg, 


force component at right angles to axis x. 


M 


kg m, 


bending moment. 


<ld 


kg nf x , 


outside force at right angles to bar axis 


q l 


kg m~ \ 


outside force parallel to bar axis. 


m 


kg, 


moment loading of bar elements. 


r 


m, 


radius of curvature of the elastic line. 


a, 




slope of the elastic line. 


X 


m, 


deviation of bar axis from straight line in 
unloaded condition. 


€ 


m, 


"amplitude" of bar axis in unloaded condi- 
t ion. 




m, 


deflection. (See fig. 2.) 


£ 


m, 


shifting in direction of x. 


S=/ P 


dt kg s , 


shock impulse, (See footnote, page 26.) 


E 


kg m~ 2 , 


elasticity modulus. 


I 


m, 


length of bar* 


i 


m, 


radius of gyration of bar section. 
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P m 3 , area of tar section, 

J in 4 , inertia moment of "bar section. 

p kg s s ia~ 4 , density. 

tf x kg nT 2 , normal stress in x direction. 

€x»..- elongation in x direction, 

t . s, .time interval. 

T s, period of oscillation of the free funda- 
mental oscillation. 

T s , , shock period, 

p s~ l , frequency of oscillation. 

v - m s" 1 , velocity of sound in "bar material. 

Cp , . phase lag, 

a, ratio of shock load to Eulerian load, 

b, , ratio of shock period to oscillation period 

of the free, transverse fundamental oscilla* 
t ion, 

c, ratio of the maximum moments (taken abso- 
lute) in the static and dynamic case. 

\, proper values. 

A,B t C,D f k x f k 3 , constants. 

Indices : 

n=l,2,3, the natural numerals. 

0 refers to quantities appearing with ten- 

sion - 0. 

E, quantities representative of the Eulerian 

buckling. load. 

., .., etc. , dots over a symbol denote its 1st, 2d, etc. 

derivation in time rates, 

~, ^, etc., dashes over a symbol denote its 1st, 2d, 

etc. derivation with respect to a length 
(x). 
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III. THE DIJfSREHIIAL EQUATIONS OP THE SYSTEM 

Within the curvilinear system of coordinates x y 
the equilibrium equations for an element of the "bent "bar 
of length ds are as follows (fig. l): 

( L + ds^ cos dp- (v +'|B- ds") sin dcp-L+q; ds=0 U) 
\ cs / \ ds / 

+ || ds) sin d£p+ + || ds) cos dcp-D+q d ds=0 (2) 

( M + H ds ) + ( L + ds ) ds sin l $ + 

+ (D + || ds) ds cos ^ - M - | qi ds 3 sin ^ + 

+ | q d ds 2 cos f + in ds = 0 (3) 

With r d cp = ds, wherehy r = curvature radius, 
these equations, upon dcp -~->0 and disregarding the infi- 
nitely small quantities of the 2d and 3d order and with 

sin dcp ~ dcp 

cos dcp ~ 1 

reduce to 

|i. liJt^O (la) 
|L + ||+q d = 0 (2a) 
+ D + m = 0 (3a) 

ds 

The deviations of the l)ar axis from a straight line 
in unstressed condition as well as the deflections in the 
processes analyzed hereinafter, are assumed small compared 
to the length of the tar, and the choice is a rectangular 
system of coordinates x y such that in first approxima- 
tion axis x coincides with the axis of the har in un~ 
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stressed attitude. Then the relations between the force 
components in the curvilinear system x y and the force 
components in the rectangular system x,y are 

L = X cos a + I sin a 
D=~Xsina+Ycosa 
or, since a « sin a = tan a - y 1 , 

4a = i- =: v» 
dx r 7 

| 5 = cos a = 1 
ds 



As a result: 





I = X + 


Y y' 






D = - X 


y« + 


Y. 


3L = 

ds - : 


g + 3i 


y« + 


Y y" 


3D _ 
ds " 




, 3Y 

+ 


- X y» 


8M _ 


oli 






are 


wr itt en 


into 


(la) to (3a), whereby, 



librium equations for the slightly bent bar become: 

- || y« + || + q d = o (2d) 

|| - X y' + Y + m = 0 (31)) 
By eliminating Y they reduce to 
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|$- q* 7« + = 0 (4) 



3x 

Assuming zero outside load at the "bar element, i.e., 
q\ and q^ to he mass forces of the bar element and m 
the mass moment of the "bar element q^ y 1 are negligible 
relative to qi , "because, first, in sufficiently thin 
bars the oscillation frequency and through it the mass ac- 
celeration in transverse direction is small relative to 
the corresponding quantities in the longitudinal direc- 
tion; second, y* is a small quantity according to the 
premises. 

Thus the differential equations read: 

§| + q t = 0 (4a) 

SIB. x y « . , 4+ g. 0 (5a) 

Now y is the deviation of the bar axis from straight 
line in unstressed condition and r\ is the deflection, so 
that y + T| must be substituted for y in (5a) (fig. 2). 

X and M are expressed in terms of deformation: 
X=a x F=EF £j = E ? g~ 

M = E J 

and the mass forces and q^, and the mass moment m 

as b 



q % = - P J 



at 3 
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Herewith (4a) and (5a) "become:. 




■ o 



(4b) 



(5b) 



t,j • These two equations (4b) and (5b) constitute the dif- 
ferent ial •' equat ions of the buckling stress due to shock 
"of a bar with constant cross-sectional dimensions. 



IV, SOLUTION FOR A BAR HINGED AT EITHER END 

f 

1. The Eree Longitudinal Oscillations of the Bar 



We repeat the well-known formulas for the free longi- 
tudinal oscillations of a straight bar. They are obtained 
by resolving (4b) conformably to the generalized equation: 

I = (A' sin \x + B cos \x) (k x sin p 0 t + k 2 cos p 0 t) 

The bar is assumed to be clamped or fixed at one end 
but left free to move longitudinally at the other. Then 
the boundary equations are: 



3=0 because of the first boundary condition, thus 
reducing the equation to 

1 = sin \x (ki sin p 0 t + k 3 cos p 0 t). 

Double differentiation according to x and t affords 



e 

li 

3x 



X 

EE 



0 » 



for 



x = 0 



x = I. 




which, written in (4b) gives 
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- X 2 35 + p Q 2 p = 0 p Q = X A / r y= V X, 

where v = y ™ is the velocity of sound in the material. 

Furthermore, because of 

= X cos \ x (ki sin p t + k 2 cos p t) 
together with the second "boundary equation we have: 
X JL 3tl 2n - JL _. 

Consequently, the frequencies are: 

- XL ,r 3tt , y 2n - 1 - „ 

P ™ " 21 V> 21 Vi V 

and the period is 

^"p "S'sv 2n~-T v # . 
2. The Free Transverse Oscillations of the Bar 
Disregarding the rotatory inertia for the case of 
X = 0, equation (5b) becomes: 

EJ S$ + p r 3 T= " 0 ( . 6 > 

The generalized solution is: 

r\ = 2 sin (p on t + Cp n ) (A n sin \ n x + B n cos \ n x + 
+ C n sh ^ x + D a ch ^ x) . 
The insertion of 

a X 4T - A n Tin and g-^r - - Pon 2 ^n 



in (6) gives 



/ E J 

Pon * V/ p"? . (?) 
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The boundary conditions for the "bar hinged at both 
ends read: 

T) = ol 

I for x = 0 and x = I, 

thus modifying (6) to 

<H = 2 -n_ - E A sin (p t + <p n ) sin \ x (7a) 
n n 

whereby 

X _ TT 2TT OTT 

» - r» t T 

By virtue of (7) the frequencies are: 

tt 2 Aij 4tt 3 AeT n 3 ^ /liJ _ n 2 TT 2 v , 

p on - ~fJ pF' ?v pp T"V " 'T~ v 

where i = radius of inertia of the cross-sectional area. 
As a result the oscillation period of the (n-l)th higher 
oscillation is: 

_ 3 TT _ 5j!_ 

011 ~ Pon ~ ? ~TT v i * 

The ratio of oscillation period of the transverse and 
longitudinal oscillations for the fundamental harmonic 
(n = l) in the "bar hinged at both ends to that of the "bar 
left free to shift longitudinally at one end, is: 

JLiL- JL ^ JL 1 . 

tt v i 4 1 2tt i " 

For slenderness ratios j > 2tt the transverse oscillations 

are consequently of lower frequency than the longitudinal 
oscillations. 

3. Transverse Oscillations Due to Constant Shock Load 

Eere we analyze the specific case of a "bar hinged at 
both ends being stressed under constant shock load X dur- 
ing a shock period t # We idealize the case "by disregarding 
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the longitudinal oscillation, i.e., assume the tension 
(normal force) to "bo identical at every point of the "bar 
during the period of shock. With the neglected rotatory- 
inertia the resolvable partial differential equation is, 
according to ( 5h) : 

The deviations of the "bar axis from a straight line, 
i*e., the original "bar form, are expanded in Fourier se- 
ries conformably to the local function sin X n x of (7a) 

y = Z sin X^ x = €i sin X x + £s sin 2 X x + .... + 

+ € n sin n X x + (9) 

where X » ^ is the first proper value of the free trans- 
verse oscillation (see section IV, 2) and e n = constant. 
Equation (8) resolves to 

r\ = Z ri tn 'Hxn - ^ti sin X i x + 'Hta sin X 2 * + + 
n 

+ r) tn sin >^ x * . . 

- T|tx sili X x + sin 2 X x + «... + 

+ r^ tn sin n X x + . . . . , 
where r|t = f(t) 

Putting these values of y and T) in (8) gives 
E J X 4 2 n* r| tn sin n X x + X X 2 Z n 2 r\ tn sin n X x + 



d 



+ p F Z S in n X x + X X 2 Z n 2 c n sin n X x = 0. 

a dt 2 n 

The equation for the fundamental equation is: 

E J \ 4 T) tl + X A 2 r| tl + p F HJi. + X X €i = 0 (10) 

. — _ — dt 
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and for the (n-l)th higher, oscillation: 

d 2 n 

n 4 EJ* 4 -ntn^ 3 * \ a mn+P* + ^ x ^ C n=° 

dt 

These equations are V&lid during, the actual shock pe- 
riod. After the shock the bar executes free oscillations 
whose initial conditions are contingent upon the deflec- 
tion and rate at the termination of the shock. 

a) The fundamental oscillation 

The Sulerian buckling load for the bar hinged at both 
ends and subjected to static load (see also V) is: 

P E =EJ^|-=EJX 3 (12) 

If 

By denoting the ratio of shock force to Eiilerian load 
with a, or in other words, presume 

a = - I, (13) 



,X = - a.E J a 2 and (10) becomes 
p p SLJQt + B j x 4 (i - a ) <n t = a € E J X (14) 



^ 2 

= — a . E J 

i 2 
dt 

index 1 being omitted for simplicity. 

The resolution of this equation must differentiate 
between three cases: 



1. 


X > - 


P E , 


that 


is, 


a 


< 


1 


2. 


X = - 


P E , 


n 


ti 

• 


a 




1 


3. 


X < - 




ii 


ii 


a 


> 


1 



ft) Shock load lower than the Eulerian load* 

(a < 1) . 



♦This also, includes all cases with negative a, that is, 
the cases of shock stresses in tension.' 
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Putting r\ t = f\ t + r\ t in (14), whero *H denotes th 
resolution of the homogeneous equation and r\ t the effect 
of the disturbing function, 'the resolution of the homogene- 
ous equation 2 

d T\ 

p F — + E J X 4 (1 - a) <n t = 0 

gives f[t - ki sin p t + k 2 cos p t, 

-because -- X* (1 - a) > 0. Consequently, 

pF 

p „ t 

The insertion of these values in the homogeneous equa- 
tion gives the frequency. 

p ■ * ffi^ 1 - *° ^ r -~ (X5) 

where p 0 = frequency of the free fundamental oscillation 
(see IV, 2). 

The effect of the disturbing function on the right- 
hand side of (14) is found from 

E J X 4 (1 - a) th = a € E J ^ 

at 

"» ■ r?T £ - 

Consequently , 

a . * 

r\ t = k x sin p t + k 2 cos p t + J - -'" e ^16) 

whence, 

r\ t - p (kj cos p t - k 2 sin p t) (IV) 

Assuming Tit = 0 and r\ t ~ 0 for the start of the 
shock, t = 0, (16) and (17) give 
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The solution is: 

T1t = G (1 - COS p t) (18) 

Jm Q, 

and consequently, 

^ = r—— e p sin p t (19) 
x — a 

Lot T = shock period. Then at the end of the shock: 
T\j = — ~ — e (1 - cos p T ) (20) 

T) T = 7 ~~ — e p sin p t (21) 
1 1 ~ a 

The incipient free oscillations following the end of 
the shock are, according to (6): 

r\ = T| t ri x = C sin (p 0 t + cp) sin \ x (22) 

with <P = phase shifting. For this period of the processes, 
(20) and (21) are the initial equations. Thus, 

C sin (p 0 T + q>) = € (1 - cos p T) 

x — a 

a 

C p 0 cos (p 0 t + cp) « i~"a C P Sin P T * 
The addition of the squares of these equations gives: 

C 2 - 7T-™-T*r c 3 [(1 - cos p t) 2 + -?1 sin 3 p T] 
(1 - a) p 0 ^ 

or, with due regard to (15), 

c - & — € /"i - 2 cos p T - a sln^p t. 
1 - a v 

We use the oscillation period of the free oscillation 

2rr 

T = ~ as time scale and introduce 
Po 

^ m 2tt ^ 
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so that p T zz 2 rr h ,/l~~~ a, and 



C = — * — £a/2~2 cos(3 tt h^l-a) - a sic? (2 rr W 1-a) . 
1 p * a 

Our interest centers about the maximum "bending moment 
whiqh is proportional to the maximum curvature (3 T)/dx 2 ) . 
Formula (22) concedes 

2JD. = - c \ sin (p 0 t + Cp) sin \ x. 
dx 3 



The curvature is maximum in the center of the "bar 

TT 3 



(x = \) and amounts to 



raax |4L = C X 3 = G Eg. - € ^ f (a,!.) 
B * I I 

where 



f(a,t>) = ~ J2-2 cos(2tt!) A/l-a)-a sin 3 (2ttT> (23) 
x** a 

This value is decisive for the maximum moment after 
the actual shock period. However, it may happen that a 
still greater moment occurs during the actual shock period. 
Such is tho case vhen the shock lasts at least long enough - 
until the highest possible deflection has been reached once. 
According to (18), the greatest possible deflection during 
the shock occurs once when cos pt = -1. Consequently, if 

T = H zz JL 1 
> P Po fl~-~l[ 

_ _ TT 1 1 



> PT fT'- a 2/1 - a 
the maximum moment is already reached during the shock. 
On the other hand, according to (18), 

^ = — e (1 - cos p t) sin ^ x 
1 - a 



during the shock f hence 
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^ r c\ a 3 'i 

max - 2 = 2 : c K = e t f (a, b) , 



is approximately proportional to the maximum moment, whereby 

f(a,b) = 2 — (24) 
1 - a 

* V ....... ........ 

To sum up: Tor shock periods b < — j=rrrr =s. 9 formula (23) 

2 a/ 1 — a 

is valid; for shock periods b ^ — —1— =u (24) is the de- 

2J1 - a 

cisive quantity for the maximum moment. For b = — y=r^=-r^z, 

2/1 - a 

(23) steadily resolves to (24). The latter represents an 
absolute maximum value of f(a,b), which may not be ex- 
ceeded with any shock period, 

(3) The shock load equals the Eulorian load (a = l). 

In this case (14) reduces to 

£nt . !£x 4 . c . 

dt pF ... 
Integrating twice gives: 

TPT 4 S 

r\t = ~rz: \ it + ki t + ka. • 

The initial conditions: 

r\t = 0 and r\ t = 0 for t - 0 
concede k a = 0 and k 2 = 0 • Therefore, 

and for the end of the shock, 

and , VT 4 

Tfc = H X € T (27) 
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The time interval after the actual shock is again com- 
puted with (22) and the initial conditions for this shock 
period are posed in (26) and (27). Thus, 

C sin (p 0 T + <p) » JSJ-V c T s 

C p 0 COS ( Po T + <p) = M \* C T 



hence . 



0 = — K £ T 



jl. 

V 4 p n 3 



Po 

or, when taking 

3 



t = §31 "b and t = p 0 



into consideration. 



C = 2 rr b c ^/ tt 2 b 2 + 1. 
Equation (22) again yields 

max 0 = C X 2 = € Hi f( a ,b), 

where ^ 

f(a,b) = 2 tt b fV* b 2 + 1 (28) 

Contrary to the case of a < 1, equation (25), appli- 
cable during the shock, is now aperiodic, hence the deflec- 
tion during the shock may not exceed that at the end of the 
shock. But the latter is, conformable to (28), the start 
of the free oscillation. 

7) The shock force exceeds the Eulerian load (a > l) 

In this case, 

« X* (1 - a) < 0, 

and the solution of the homogeneous equation corresponding 
to (14), manifests: 
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f\% = k x sh p t + k 3 ch p t f 
whence, 2 ^ 

which, written in the homogeneous equation, gives 

p 2 = g£ \ 4 (a - l) = p 0 3 (a - 1) (29) 

p 0 = frequency of the free fundamental oscillation. 
E J X* (1 - a) fj t = a C E J ^ 



concedes the particular integral at 

a 1 



Consequently t 

rit = k i sh p t + k 2 ch p t - a~~T € (30 ^ 
T| t - p (k x ch p t + k 2 sh p t) (31) 

Assume Tfc = 0 and T| t = 0 for t ~ 0, Then (30) 
and (31) give 

a 

k 3 = r- C and k x = 0 

a x 

whence, 

•Ht = — * e (ch p t - 1) (32) 
a l 

and for the end of the shock t = T: 

^ T = a"- T C (ch P T " 1} (33) 
TJt = a~-T ^ P sh p T .' . (34) 
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For the free oscillations after the shock (22) and 
the initial conditions (33) and (34) are again applicable, 
ffhen, 

C sin (po T + cp) = — - € (ch p t - 1) 

a - 1 

C Po COS (p 0 T + Cp) = — ™- € p sh p T 

and, with due allowance for (29) 

0 = — e J~2 - 2 ch p T + a sh 2 pT 



or with p P 

p T = 2 TT "b = 2 TT *b a/ a — 1, 

Po 




According to (2.2) the maximum curvature is again 

max ft = C X 2 = € % f( a ,b) 

ox i 

where 

f(a,b) = ~ ch(2TTb^~l) -ha sh 2 (2rr"b/a-l) (35) 

Since (32), applicable during the shock, is aperiodic 
for r)t and increases with t, no greater deflection can 
occur during the shock than the maximum deflection reached 
after the shock; but that is comprised in (35). 

b) The (n-l)th higher oscillation 

Substituting X = - a E J X, equation (11) becomes: 

p j !LH|n + n 2 E j \ 4 ( n 2 - a ) Titn ^ n 2 a E J X* e n (36) 
dt 

Again we differentiate between; 

1, a < n 5 

2, a 5= n 2 

3, a > n 3 
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of which the first is of primary interest* Even with the 
first higher oscillation (n = 2), cases 2 and 3 refer to 
shock loads at least 2 =4 times as high as the Eulerian 
load. Such excesses of the Eulerian load may he disregard- 
ed and the analysis confined to a < n 2 . 

The resolution of (36) for a < n 2 is similar to that 
of (14) for a < 1« Let the solution of the homogeneous 
equation he 

'^tn = k i 3in P n 1 + k s cos Pn * 

so that g 

dt^ Pn ltn 
which gives the frequency •• 

Pn 85 * ^ /n 2 - a = n p 01 /5 r TT (37) 



. - frequency of free fundamental oscillation. (See 
section IV ,2.) 

The effect of the disturhing function follows from 

n 2 E J X 4 (n 2 - a) rjtn = a 2 a E J £ e n 

at 

^tn = ~ 2 ~~~ 
whence the solution of (8) at 



a 

a 



r) tn = k x sin p t + k 2 cos p t + s ^ B 



a 

n~ - a 



k x and k 3 are again defined from the initial conditions: 

^tn = 0 and ^tn = 0 for t 53 °" The resttlt 
a 

"Htn = "a c n (1 - c °s P n *) (38) 

n" - a 

"Htn = IT P » G n sin Pn* = TT'~ p oi € » sin Pn* 

- a V n"- - a 
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For t ■« t, we have : 

a 

ri Tn = — € n (1 - cos p n t) (39) 

n ~ a 

• Ha 

'Hth = - y — Vox € n sin Pn T < 4 °) 
V n 2 - a 

After the shock the "bar executes the (n~l)th free 
higher oscillation, which is governed "by 

r\ n « 0 n sin (n 2 p 01 t + Cp n ) sin n \ x (41) 

With C and cp defined from (39) and (40) as initial con- 
ditions: 

C n sin (n 2 p 01 T + cp n ) = '-g-? c n (l - cos p n t) 

n - a 

c n n2 Poi cos(n 2 p 01 7 + cp n ) « ~=^A-==?. p 01 c n sin p n T 

V n^ - a 

gives C at 



a f~ a 

c n = ^~r~a €n v 2 ~ 2 008 Pn T ~ n 2 sln Pn T ' 

Because of (37) it is 

p n T = 2 n it h J *- a 

when T = J-H t- Therefore, 
Poi 

C n = — c n / 2-3 cos(2nTT"b\/ S n 2 -a)- 4r sin 2 (2nrrh v^li^-a) . 
n -a v n 



According to (41) the curvature is 



|JD§ = _ Cn n 2 x 2 sin (n 2 p ox t + cp n ) sin n X x 



and the maximum curvature. 



- 2 2 



max g~^r - C n n s A, = -3- £ n f(n,a,b), 
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where f(n,a,*b) = 

= fz-Z cos(2nTTb/n 2 -a) - ™ sin 3 (2nnb Jxf^a.) (42) 

n ~a n 

The validity of (42) for the maximum moment is decis- 
ive only when it occurs after the shock* If the actual 
shock lasts at least long enough to permit once the occur- 
rence of the highest possible deflection, then the maximum, 
moment' occurs during the actual shock. This is the case 
according to (38) whon cos p n t = - 1 , that is, when the 
duration of the shock is 

T g JL = -JL -J 

Pn 11 p oi /n 2 - a 

2 n a/ n 2 - a 
In this case the curvature is 

(3 'n a s 

_ — § = -g- % n 2 >v (1 - cos p n t) sin n X x 

ox n - a 



according to 

a 

n^ ~ a 



T| n = ™ €(i-cosp n t) sin n X x 



and the maximum curvature is 

max I = _^ e a n s \ = €n f( n ,a,l>), 

o x n - a i e 

with 2 a n 3 

f(n,a,b) = (43) 
n - a 

Summed up: for shock period 

D <. ~ — =ssr r 

2 nvn 3 - a 

equation (42) is applicable; . . . 
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f 6r 

^ > 1 

2 a 3 - a 



(43) is valid as the quantity deciding the maximum moment. 

1 

2 nvT 2 - a' 

ter represents the maximum value of f(n,a f b), which is 
not exceeded in any shock period. 



When b = "" 7"? ~ > ( 42) "becomes (43). The lat- 



V, NUMERICAL INTERPRETATION AND DEDUCTIONS 

The "behavior of the bar during shock load is best eval- 
uated by comparing it with its behavior under static load. 

The differential equation of the static load is: 

E J J-J - X r\ = X y (44) 

with y given from (9)« Limited to the first term of (9), 
the resolution gives: 

r\ - C x sin X x + C2 cos X x 

or r| = C x sin X x (45) 

since r\ = 0 when x = 0. 

Putting (45) in (44) gives X= ^ and 

Ci = — i± — ■ ■ ■ . 

E J TT 2 , 
-3 1 



consequently, 



X I 



sin X- x 
Tj = 

X 1 



(46) 



for - X = P E - E J , 



24 



IT.A.C.A, Technical Memorandum ITo. 748 



that is, the Eulerian load, r\ = «>• Therefore, and with 
consideration of (13), 

T\ = — - — € x sin \ x 
1 - a 

whence a 3 ^ ^2 

max ~^ ^ C X f( a ) (47) 

with f(a) = — — ' 
1 - a 

The ratio c of the moments due to shock load and 
static load is: 



c - J 2 - 2 cos(2iTbyr - a) - a sin 2 (2rrb Jl - a) (48) 

according to (23) and (47) when a < 1 and the duration of 
the shock is 



Z jl - a 

and c = 2 (49) 

according to (24) and (47) when the shock period is 

1 



b ^ 



> 



2</l - 



a 



It is readily seen that the c terms are dependent on 
the magnitude of the eccentricity c. (See fig, 3.) 

It will be noted that the ratio c of the dynamic and 
static stress for short shock periods is smaller, at longer 
periods greater than 1 and its maximum value 2. Moreover, 
for equal shock dv.rations, ratio c is smaller as the 
shock load is higher. 



When plotting, as in figure 4, that shock load b 
against shock load a for which the dynamic equals the 
static stress, it is seen that comparatively long shocks 
are necessary vicinal to the Eulerian load to raise the dy- 
namic stress on a level with the static stress* With a 
shock equivalent to 0,97 times (approximately) tho Euler- 
ian load, this shock period equals tho natural oscillation 
period of the froe bar. 
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If the form of the bar is such as to oxactly produce 
: .the (n-l)th higher harmonic, by the same argument the nth 
stability limit is 

- X = EJ ' 

when the "bar produces only the fundamental oscillation and 
.•the decisive function for the maximum moment is 

f(a,n) = -Si-*- (50) 
n 2 - a 

When 

2 n</n 2 - a 
the comparison with (42) and (43) gives 

c = J 2-2 cos(2nTTbv / n 2 -a) - -% sin 2 (2nTTb ^/^-a) (51) 

n 

and e = 2 

when 

D j=z — ^r=L • 

2 n fj n 2 - a 

The range of validity of (51) is limited to very short 
shock periods. Even for n = 2 and a - 1 , the upper lim- 
it of "b is 0.144 only. Since for very short shock peri- 
ods the premises of the calculations are in any case hard- 
ly met (reference l) , the evaluation of (51) may he fore- 
gone, especially since the case where exactly only the 
(n-l)th higher oscillation occurs, is practically without 
s ignif icance. 

By contrast, the case where the "bar shape is such as 
to incur several oscillations concurrently, is much more 
important. But obviously this case does not lend itself 
to general treatment, "because the results are substantial- 
ly affected by the relative magnitude and the sign of € n 
in (9). 

"The resolution of (44) is applicable only to the cases 
for the evaluation of the data of the dynamic investiga- 
tion in which the load lies below the stability limit, i.e., 
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for n = 1 in the a < 1 range; for n = 2 in the a < 4 
range, etc. For loads above the stability limit (44) should 
be disregarded in favor of the more exact equation of the 
elastic line, the resolution of which is, however, quite 
complicated, For that reason the comparison of the data 
is limited to n = 1 for a < 1, a = 1, and a > 1, re- 
spectively. The results are illustrated in figure 5, in 
logarithmic scale* The functions f(a,b) from (23), (24), 
(38) /'and (35), proportional to the maximum moments are 
plotted for divers a against b» 

It is readily seen that the Eulerian load may be ex- 
ceeded in buckling stresses due to shock, provided the 
shoe?*: period itself is short enough, Koreover, the maxi- 
mum dynamic stresses are fostered by increasing shock load 
a and period b. But, while attaining a limit value for" 
shock loads below the Eulerian load (a < l) for a given 
duration of shock, which cannot be excoedod in any shock 
period, they increase arbitrarily at shock loads above the 
Eulerian load. 

VI. EFFECT OF SHOCK IMPULSE / F dt ON THE STRESS 



Frequently it is assumed that the stress due to shock 
load p is dependent only on the shock impulse / P dt , 
that is, individually unaffected by the magnitude of the 
shock load and the duration.* 

The results of the present paper disclose the error 
of this assumption, for otherwise only the product a b 
would appear as sole variable of the terms for f(a,b). 
Again, it ma;' be asked whether or not it would be approxi- 
mately correct. For that reason, we compute the functions 
f(a t b) versus b for several values of S = / P dt = ab. 



*'It is common practice to designate the time integral 
/ P dt as "shock load," whereas the quantity P is not 
specifically expressed. This practice is probably due to 
the concept that only the time integral / P dt is decis- 
ive for the shock process, whereas no special importance 
attaches to quantity P. But the author s of this paper 
have, on the strength of thoir investigations, drawn dif- 
ferent conclusions, and believe it, in fact, to be more 
logical to express P, which has the dimension of a force, 
as "shock load" and time integral / P dt with the dimen- 
sion of a force times time interval as "" shock impulse." 
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Concentrating on the fundamental oscillation, there 
are three ranges of "b for a stated value of S = ah: 

li range 0 < b < b x for a > 1. 

2, » b x < b < h 3 for a < 1 to the 
extent that the maximum stress occurs 
^.?XS3L the actual shock period. 

3» range b 2 < b when the stress occurs 
during the actual shock period. 

The values for b x and b 2 may he defined as fol- 
lows: 

a = 1 for h lf consequently, h x = S; 

*s = —7== for b 2 , 
2 V 1 - a 



or 



"because 



t 2 = I (s +ys^ + 1) 



s 

a = ~. 

o 2 



Tor f(a f b): 

equation (35) is valid in the range of 0 < b < hi 

» (23) » » " " » " bi < h < b 2 

ii ( 24 ) ii fi ii i. i» ii bg < t 

ii ( 28 ) ii u ii ii n ii _ 

The value of f(a f h) for b = 0 is obtained by put- 
ting a - S- in (35) and permitting b to approach zero* 
Then, 



lim f(a,b)= / ! sh 2 (2 tt a/s b - b 2 ) = 
b-0 ^ b 



, „ s ^^^L^J^^JL=^ = 2 n S. 

S b - b 2 
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With these formulas we computed f(a,b) for S = 0.25* 
0.50, 0.75, 1.00, and 1.25, as well as for various b val- 
ues. The results are plotted in figure 6. The ordinates 
of the individual curves are noticeably not approximately- 
constant with the parameter S or, in other words, the 
stress due to a shock load cannot even he approximately 
given in function of the shock impulse. On the contrary, 
shock load and shock period must he individually known 
if the shock stress is to he determined. 

This result is of great importance for shock tests. 
Shock "load and shock period must "be included in such ex- 
periments, although this will ho more difficult to accom- 
plish than recording the shock impulse. 



APPENDIX 

Effect of Minor Changes of the Original Bar Shape 

on the Results 



The premise of the interpretation of the results was 
the selection of the original "bar shape such as to precise- 
ly insure the occurrence of the fundamental oscillation 
due to shock load. In the example hereinafter, we attempt 
to show the effect of a minor chango in the original shape 
of the har on the results. For simplicity the range is re- 
stricted to a < 1, that is, to the range within which the 
ratio c of the dynamic and static stress is readily ob- 
tainable. 

We assume the shape of the bar such as to develop 
aside from the fundamental oscillation, yet the second 
higher oscillation. Then (9) reads: 

y = € x sin ^ x + € 3 sin 3 ? x 

with — ~ \ 9 

Tor shock periods, 




the equation 

T, > * = 1 „ 

2n/n 2 - a 6 *J~9 - a 
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is particularly applicable. According to IV, 3a and IV,3"b 
the maximum deflections occur during the shock. The cur- 
vature of the elastic line is expressed "by 



9 e 3 
+ 9--T 



- — - x - (l-cospt) sin X x + 
1 - a 

(1 - cos p 3 t) sin 3 X x J J 



(52) 



The points of the maximum and minimum values of this 
function result from the resolution of 



15a 

dx 3 



a X° 
+ 9 ""a 



— (l-cospt) cos A x + 
x a 

(1 - cos p 3 t) cos 3 X x J I = 0 



(53) 



With 



cos 3 X x = 4 cos 3 X x - 3 cos X x 
(53) resolves to cos X x = 0 

and 



(54a) 



*- (l •-cos p 3 t)(4 cos X x— 3)=0 ( 54"b) 



e , / , 27 € 

11- cos pt) + 

1-a 9-a 

First we consider the maximum for x = of (54a). 



Here (52) gives: 



a r\ 



a X 



9 £ 5 



; (l-COS pt) - (l-COS p 3 t) 

1-a 9-a 

9 e 3 

To the extent that |e 3 | = €j, hence 



(55) 



1-a 



9-a 

the limits of the maximum curvature are 



2a xV « - |JE^ 

\l-a 9-a/ 



9 €. 



3 2 jq 
Sx" 2 



max 
x= 



2a X 2 - C JL 
1-a 



(56a) 



when € x and e 3 have the same sign, and 
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2aX 2 i- 1 - 
1-a 



.J3 

ax 3 



max 

x=4 



Vl-a 9- a / 



( 56b). 



when c x and c 3 have different signs. 

• Comparing these results with the value of the curva- 
ture at x = 1/2 in the static case 



-2 

d r\ 




"rran 
dx 





^(rK--r^)|- 



the ratio c of the dynamic and static stress is found to 
lie "between 

t x (9 - a) 



2 ^ c = 2 



€i (9 - a) - 9 € 3 (1 - a) 



(57a) 



c ^ 2 



when c x and e 3 have the same sign, and 

€ x (9 - a) - 9 e 3 (1 - a) 

when they have a different sign. 

This leaves the question, whether or not at some 



(57b) 



may 



point other than x = g an upper limit of 
occur which exceeds the values given in (57a) and ('57"b) . 



Bin 
a? 



max 



Superior limit values of 



ax 2 



max 



may occur at 



points other than x - * only when x meets equation 

(54h). Let x x "be a real root of this equation. Then 

P \ 3 1 9 - a 1 - cos p t 

| A. *»*- = - — *- 



COS' 



*1 



hence , 



4 108 1 - a €3 1 ~ cos p 3 1 
sin 3 X x 1 = sin X x x (-l+4 cos 2 X x x ) = 

= sin X x,(2 - - 1 - li. ilSaLJLl) 

x ^ 27 1-a £ 3 1-cos p t/ 

Then, according to (52) the curvature at x = Xi is: 







3^ 
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a X sin X x x 



2 _ 

3 1 "- a 



+ _ ?_ d _ cos p 3 t)J 



9 - a 

The upper limits of (58) are: 



(58) 





< 




x=x x 



2a X sin X x x ( - -— x — + - — — ) 
x \ 3 1 - a 9 - a/ | 

when £i and e 3 show the same sign, and 

4 

3 



(59a) 



^ 2 


< 


Bx 3 


X=Xi 



- A. *- S xn A, x, 

1 - a 1 



(59b) 



when otherwise. 



The requirement that the curvature at x = k 1 for 
equal signs e x and € 3 shall at the most be the same 
as the maximum value of (56a) results in 



2a X 2 sin X x x (% + ^) 

\3 1 - a 9 - a/ 

This condition is always met: 



2a X 



1 - a 



i.e., for 



because 



2 e x 18 c 3 

3 1 - a 9 - a 



< 1 
6 



1 - a 



(60) 



1 - a < 1 
9 - a = 9 ' 

For different signs of e x and c 3 the result is 
similar, according to (56b) and (59b): 



3 X r~a Sin X * l 



2a X 3 (-Ji-.-LiLM. 
\1 - a 9 - a/ 



This condition is met for every value of € 3 as can 
readily "be seen. 
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Figure 3. -Ratio c of maximum moments under dynamic and 

static load versus ratio "b of shock period 
to period of free oscillations for various ratios a of 
longitudinal force to Eulerian load. 




Figure 4. -" Shock period" "b for divers "shock loads" a, 

for which the maximum dynamic and static 
stress are equal, i.e. c = 1 . 
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Figs. 5,6 



1.05 



a=,75 




EEliOO Eulerian load 



0.1 0. 



ZTigure 5.- Quantity f(a,h) 
proportional to 
the maximum moments plotted 
against "shock period" h, 
for various "shock loads" a. 
(logarithmic scale) 




